On the detection of condensate vortex states 
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We discuss a detection scheme that reveals the existence of vortex states in a cylindrically sym- 
metric condensate trap. It relies on the measurement of the second-order correlation function of the 
Schrodinger field, and yields directly the topological charge of the vortex state. 
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The rapid progress in the experimental generation and 
manipulation of Bose-Einstein condensates in low density 
trapped alkali vapors fl|-|5| opens up the way to the de- 
tailed study of the thermodynamic and dynamical prop- 
erties of weakly interacting quantum-degenerate gases. A 
topic of much current interest is the superfluidity of these 
samples. Superfluidity is related to the rotational prop- 
erties, and in particular to the existence of vortex states 
in quantum gases ||. While no vortex states have been 
launched in low density atomic condensates so far, their 
existence has been numerically studied 0-^) and their 
stability analyzed {f(J] . Optical methods to launch vortex 
states were recently proposed in Ref. [p 12|, and Jackson 
et. al. jl3| have produced numerical solutions of vortex 
formation in a weakly-interacting condensate by piercing 
it and subsequently slicing it with a blue-detuned laser 
at a velocity exceeding the critical velocity. 

In view of these developments, it is therefore timely 
and important to discuss ways to detect vortices in 
trapped atomic condensates. One possibility suggested 
in Ref. Jll| is based on measuring the spatial absorption 
profile of the sample. This method is sensitive to the ro- 
tation of the condensate via the rotational Doppler shift 
jl4| , but does not give a direct measure of its topological 
charge. Rather, the quantization of the vortex motion 
must be inferred from other condensate parameters such 
as number of particles, system size, etc. In contrast, the 
scheme that we propose and analyze in this note produces 
a spatial interference pattern that results directly from 
the nonzero atomic angular momentum of the vortex and 
gives a direct measurement of its topological charge. In 
addition, a slight variation of that method allows a direct 
demonstration of the existence of persistent currents in 
the sample. The method is based on the ionization de- 
tection scheme discussed in |lq ], whereby one measures 
normally ordered correlation functions of the Schrodinger 
field operator \P(r, t). Specifically, we show that the two- 
point second-order correlation functions provide direct 
information on the persistent currents and the topologi- 
cal charge of the vortex state. 

The proposed schemes [^1 12 to create vortices in 
cylindrically symmetric traps rely on Raman transi- 
tions between two hyperfine levels of ground state alkali 
atoms, the effective coupling resulting from virtual elec- 
tric dipole transitions induced by far-off resonant laser 



beams. For this discussion we assume that the system is 
effectively two-dimensional as a result of strong confine- 
ment in the z-direction. A coherent rotational coupler is 
realized if the lasers used for the Raman transitions carry 
angular momentum associated with their field envelope, 
which can be transfered to the center-of-mass motion of 
the atoms. This can be achieved e.g. jll| by placing the 
condensate at the focus z = of two Laguerre-Gaussian 
beams of polarizations a + and cr - propagating along the 
z direction and with Rabi frequencies 
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where p 2 — x 2 + y 2 , w is the beam spot size, and ip is the 
azimuthal angle. The transverse structure of these fields 
is characterized by the radial mode number p and the az- 
imuthal mode number m fl6| . In particular, they have a 
phase singularity at their center with topological charge 
±m and associated angular momentum ±hm. As a result 
of the Raman transitions the laser fields couple a conden- 
sate in the ground state (r\g) = ip g (r)\F = 1, Mp = —1) 
with no angular momentum, where tpg( r ) is the ground 
state wave function of the trapping potential, to the vor- 
tex state (r\v) = ip v (r)\F = l,Mp = 1) with 2hm angu- 
lar momentum per atom and topological charge 2m. For 
the geometry discussed in Ref. fTl|| and using m = 1 the 
vortex wave function is given explicitly by 



^ v {v)=Np 2 e^H g (p) 1 
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with J\f a normalization constant, and 8(r) — 2ip. The 
key point from this discussion is that the vortex state 
created in this manner has a topological charge 2m and 
associated azimuthal variation exp(2im<p) in the general 
case. 

Since (v\g) — 0, this system can conveniently 
be described in terms of the spinor <p(r, t) = 
col(<frg(r, t), <j> v (r, t)). Then in the mean-field Hartree ap- 
proximation, and assuming that many-body effects do 
not significantly alter the spatial profiles of the ground 
state and vortex state so that we may restrict ourselves 
to these two modes, the components of the spinor evolve 
according to 
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&,(r,i) = (3(t)Mr), fl (r,i) = a(iM,(r). (3) 

For the case m = 1 and f3(t) are governed by the 
equations of motion (5) and (6) of Ref. [fn], and this 
serves as a concrete example. 

More generally, in the following we consider a generic 
situation where the external perturbation responsible for 
the coupling creates a quantum degenerate gas in the su- 
perposition (El), but with the phase 9(r) of the excited 
state in Eq. (0) kept general. This will allow us to point 
out the specific signatures of a vortex state that distin- 
guish it from a non- vortex state. 

In the two-mode approximation the condensate is de- 
scribed by the two-component Schrodinger field operator 
\fr(r) = (\& fl (r), \&„(r)) satisfying the boson commutation 
relation [\E^(r), #],(r')] = S U 'S(r-r'),£ = {v,g}. Assum- 
ing then a system composed of iV-atoms, the quantum 
state of the system becomes 
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or, in terms of the mode creation and annihilation oper- 
ators 



a] = / d 3 rV^(r)*t(r), 
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d 3 r^(r)*(r), (5) 



we obtain 
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k 

(6) 

Hence, the sample is in an entangled superposition of 
the ground and vortex states, a result of the fact that 
the total number of particles is conserved, but not the 
individual particle numbers in the two states. This is 
similar to the situation of split condensates discussed in 
Rcfs. 018). 

Due to the assumed cylindrical symmetry of the trap- 
ping potential, the existence of a vortex state cannot 
be demonstrated by off-resonance imaging]^] which mea- 
sures correlation functions of the sample density p(r, t) = 
^(r, t)^>(r, t). Since this is simply the sum of the con- 
densate and vortex density profiles, and the vortex den- 
sity profile is cylindrically symmetric, the density profile 
does not reveal the phase singularity associated with the 
vortex. What is needed instead is a measurement scheme 
which involves correlation functions of ^(r, t) itself. As 



1 The situation would be different if the ground and vortex 
states corresponded to the same electronic state. 



discussed in Ref. Eq), these functions can be extracted 
in an ionization scheme whereby one or more tightly fo- 
cused lasers are used to selectively ionize atoms in small 
regions of the condensate plus vortex system. The mea- 
surement proceeds then by detecting the ionized atoms, 
which play the role of a detector field. 

We consider specifically a two-point detection scheme 
in which two ionizing laser beams are focussed at loca- 
tions ri = (pi,(fi) and r 2 = (p 2 ,<£ 2 ). For that geome- 
try, and assuming that the lasers are focussed onto spots 
small compared to the dimensions of the condensate, the 
ionization scheme measures the probability w 2 of jointly 
ionizing an atom at ri and the other at r 2 as a function of 
normally ordered Schrodinger field correlation functions 
whose explicit form is 

/t+At pt+At 
dti J dt 2 

(*Vi,*i)*V2,*2)*(r2,ii)*(ri,i 2 )> 

/t+At pt+At 
dt! J dt 2 

(^(n,^)* 1 ^,^)*^,^)*^,^)), (7) 

where ry(r) is the detector efficiency and 77(ri,r 2 ) its 
cross-efficiency, as discussed in Ref. Jig]. The first term 
in W2 is an exchange contribution, while the second term 
is the direct contribution familiar from photodetection 
theory. 

Equation (0) can be considerably simplified for mea- 
surement intervals At small compared to the characteris- 
tic evolution time of the condensate. In this case we can 
set t\ =t% = t and we have simply 

w 2 (t) = 

(A*) 2 [7 ? (r 1 ) 7? (r 2 )+r ? (r 1) r 2 ) ?7 (r 2l r 1 )]G( 2 )(r 1 ,r 2 ;i) ) (8) 
where 

G( 2 >(n,r 2 ;i) = 

($(0)|^(ri,t)*t( r2>t )*( rajt )*( ri)t )|$(o)) 

= ($(i)|*t( ri )^t( r2 )#( r2 )^( ri )|$( t )). (9) 

As follows from the definitions of the detector self- and 
cross-efficiencies rj(r) and 77(1*1, r 2 ) given in Ref. |Q the 
term in square brackets in Eq. (|8|) does not vary az- 
imuthally, which is important for the following consid- 
erations. Specializing in anticipation of our subsequent 
discussion to the case p\ = pi = P and with Eq. ([5]) this 
gives readily 

G( 2 )(r 1 ,r 2 ; i )=2A(A-l)|^(p)| 2 |^(p)| 2 |a(t)| 2 |/3(t)| 2 
x(l + cos[0(r 2 )-0(n)]). (10) 

The key feature for the present discussion is the spa- 
tial dependence contained in the explicit phase differ- 
ence [#(r 2 ) — 0(ri)]. Provided that measurements are 
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performed at a fixed time t, this dependence allows us to 
determine the existence of vortex motion, in a particu- 
larly simple way, as we show below. 

Persistent currents 

The hydrodynamic formulation of superfluidity fil| in- 
troduces the velocity v s = (h/M) V0(r), where M is the 
atomic mass and 9(r) is the phase of the superfluid com- 
ponent, in our case the vortex phase measured relative 
to the ground state phase which was tacitly taken equal 
to zero. Vortex states are characterized by the fact that 
the circulation of v s is quantized, 



6(h) - 6(r 2 ) = 2m{<px - 



(14) 



dl = 2irn(h/M), 



(11) 



where n is an integer. In order to detect the circulation 
of the vortex, it is sufficient to determine its tangential 
velocity component v v . Hence the detectors can remain 
on a circle centered on the axis of rotation of the vortex 
and we have 



™ 1 dd r - 

9 dip 
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where e v is the unit vector tangential to the radial di- 
rection. For small distances |r 2 — ri|, the relative phase 
appearing in the last term in the correlation function JTc| ) 
becomes 

0(r 2 ) - 0(ri) ~ (V0) • e v dp = (M /Ti)v v (p)dp. (13) 

For a general phase variation 6(r) the local velocity, and 
hence the current of atoms, will vary azimuthally. How- 
ever, for a vortex of topological charge 2m, we have 
v v = 2mh/Mp independent of the azimuthal position 
of the pair of closely spaced detectors, and n = 2m in 
Eq. (11). Hence, moving the pair of detectors along a 



circle while keeping their distance dp fixed allows one to 
determine the presence of persistent currents, v v ((p) — 
constant. Detection of this persistent current is a key 
signature of a vortex state. From the value of this per- 
sistent current one could infer the value of n if all other 
parameters were known. Next we describe a second mea- 
surement which yields the topological charge more di- 
rectly. 

Topological charge 

In addition, it is also possible to carry out a different 
class of measurements where detector 1 is held at a fixed 
position relative to the vortex core, while detector 2 is 
moved on a circle relative to detector 1. In that case, 
will exhibit oscillations as the relative azimuthal angle 
(ip i — </? 2 ) between the two detectors is varied, the phase 
difference being 



Thus, the second-order correlation function shows inter- 
ference fringes as the relative azimuthal angle of the de- 
tectors is varied, and the topological charge of the vor- 
tex may be extracted as twice the number of (complete) 
bright fringes in the interference pattern as detector 2 
is moved through a full circle. Note that in contrast to 
the preceding scheme, which requires knowledge of the 
system parameters, e.g. the atomic mass, to determine 
v v absolutely, the present method reveals the topological 
charge of the vortex, 2m, from a global property of the 
interference pattern, the number of bright fringes. 

To summarize, we have discussed a measurement 
scheme that permits to fully characterize vortex states 
in low density atomic condensates, yielding both direct 
evidence of persistent currents as well as a parameter- free 
determination of the topological charge. 
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